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O ; Abstract 

Q \ We calculate the maximum mass of neutron stars for three different equations of state (EOS) 
CO ' 

based on generalized Skyrme functionals that are simultaneously fitted to essentially all the 2003 
nuclear mass data (the rms deviation is 0.58 MeV in all three cases) and to one or other of three 
different equations of state of pure neutron matter, each determined by a different many-body 
calculation using realistic two- and three-body interactions but leading to significantly different 



degrees of stiffness at the high densities prevailing in neutron-star interiors. The observation of a 
neutron star with mass 1.97 rfc 0.04 Ai® eliminates the softest of our models (BSkl9), but does 
not discriminate between BSk20 and BSk21. However, nuclear-mass measurements that have 
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f\j ■ been made since our models were constructed strongly favor BSk21, our stiffest functional 
PACS numbers: 21.10.Dr, 21.65. Cd, 21.65.Mn, 26.60.Dd, 26.60.Kp, 04.40.Dg 
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The nuclear energy- density functional theory aims at describing various nuclear systems 
containing a large number of nucleons. It has been very successfully applied to the study of 
the structure and the dynamics of medium- mass and heavy nuclei lj . This approach is also 
particularly well-suited to the determination of the properties of dense nuclear matter in 
extreme astrophysical environments such as supernova cores and neutron-star interiors [2j. 
In this framework we recently constructed a family of three functionals, BSkl9, BSk20 and 
BSk21 3[, each one intended to provide a unified approach to the description not only of 
the different regions of neutron stars (including the possibility of superfluidity [4]) but also 
of supernova cores. These functionals have all been derived from Skyrme forces of the form 

Vij = t„(l + x P a )5( rij ) + + x 1 P a )^ [p% <S(r y ) + Py] 

+ t 2 (l + X 2 Pa)-^Pij-K r ij)Pij + g*s(l + X 3 P a) ™0T $( r ij) 

+ lh(l + x,P a )^ [pln(ry5(r tJ ) + 5( riJ )n(rypl] 

+ t 5 (l + xsP^-^Pij.niry 5(ry)py + -^Wofa + tr s ) -py x dir^p^ , (1) 

where r»j = r» — r 3 -, r = (r, + r 3 -)/2, p i3 - = — i/t(V, — V 3 -)/2 (this is the relative momentum), 
P CT is the two-body spin-exchange operator, and n(r) = n n (r) + n p (r) is the total local 
density, n n (r) and being the neutron and proton densities, respectively. These forces 

are generalizations of the conventional Skyrme form in that they contain terms in t± and t&, 
which are density-dependent generalizations of the usual t\ and t 2 terms, respectively Q]. 

The parameters of this form of force were determined primarily by fitting measured nu- 
clear masses, which were calculated with the Hartree-Fock-Bogoliubov (HFB) method. For 
this it was necessary to supplement the Skyrme forces with a microscopic contact pairing 
force, phenomenological Wigner terms and correction terms for the spurious collective en- 
ergy. However, in fitting the mass data we simultaneously constrained the Skyrme force to 
fit the zero-temperature equation of state (EOS), i.e., the energy per nucleon e as a function 
of the density n, of homogeneous pure neutron matter (NeuM), as determined by many-body 
calculations with realistic two- and three-nucleon forces. Actually, several such calculations 
of the EOS of NeuM have been made, and while they all agree very closely at nuclear and 
subnuclear densities, at the much higher densities that can be encountered towards the cen- 
ter of neutron stars they differ greatly in the stiffness that they predict. It is in this way 
that we arrived at the three different effective forces of this paper: BSkl9 was fitted to the 
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FIG. 1: (Color online.) Zero-temperature EOSs for neutron matter (NeuM) with functionals of 
this paper. Also shown are the realistic EOSs FP [6], WFF fl, APR Q] and LS2 [8]. Dotted curves 
denote supraluminal EOS. 

softest EOS of NeuM known to us, that of FP Kg and the one labeled "UV14 plus TNI" 

n 

in Ref. [Zfl and which we refer to as WFF, BSk21 to the stiff est, the one labeled "V18" in 

n 

Ref. 181 and which we refer to as LS2, while BSk20 was fitted to an EOS of intermediate 
stiffness, the one labeled "A18 + 5 v + UIX*" in Ref. |9| and which we refer to as APR (see 
Fig. [1]). Furthermore, the strength of the pairing force at each point in the nucleus in ques- 
tion was determined so as to exactly reproduce realistic 1 5'o pairing gaps of homogeneous 
nuclear matter of the appropriate density and charge asymmetry 10j. Finally, we imposed 
on these forces a number of supplementary realistic constraints, the most notable of which 
is the suppression of an unphysical transition to a spin-polarized configuration both at zero 



Hi- 



and finite temperatures, at densities found in neutron stars and supernova cores j^l, 5, 

The introduction of the unconventional Skyrme terms allowed us to satisfy all these 
constraints and at the same time fit the 2149 measured masses of nuclei with N and Z > 8 



given in the 2003 Atomic Mass Evaluation (AME) 12| with an rms deviation as low as 0.58 
MeV for all three models, i.e., for all three options for the high-density behavior of NeuM. 
Fig. Q] also shows the EOS in NeuM of Skyrme force SLy4 [131 ] . This force, like our own, 
was intended for use in neutron stars. We see that its EOS in NeuM lies between those of 
BSkl9 and BSk20, i.e., between the realistic EOSs of FP-WFF and APR. It was also fitted 
to a few nuclear masses, but the overall quality of the fit is rather poor, an rms deviation of 
5.1 MeV for all the measured even-even nuclei being quoted jl4l ]. 
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FIG. 2: (Color online.) Zero-temperature EOSs for neutron-star matter (N*M) with functionals 
of this paper. Dotted curves denote supraluminal EOS. 



Neutron-star structure. With each of our functionals we have calculated the zero- 
temperature EOS in each of the three regions of neutron stars. In Ref. {3] we calculated 
the EOS of neutron-star matter (N*M), the medium constituting the homogeneous core 
(assumed here to consist of just neutrons, protons, electrons and muons), which dominates 
the global properties of neutron stars. Referring to Fig. 121 we see that qualitatively, the 
EOS of N*M resembles that of NeuM for the same force, with BSkl9 still having the softest 
EOS and BSk21 the stiffest. Fig. |2] also shows the EOS of functional SLy4 in N*M. 

Aside from the obvious influence of the leptons, the main reason for the differences be- 
tween the EOS of NeuM and that of N*M lies in the presence of protons, which interact both 
with each other and with neutrons. The precision mass fit assures that these interactions are 
well represented, at least up to densities not much greater than no, the equilibrium density 
of symmetric nuclear matter (SNM). As for higher densities, the extent to which the EOS 
of N*M is not determined by that of NeuM depends on the EOS of SNM, and as can be 
seen in Fig. [3j our models are consistent with measurements of the high-density pressure of 
SNM deduced from heavy-ion collisions even though our functionals were not directly 
fitted to the EOS of SNM. Thus, all in all, we believe that the procedure for determining 
our forces has an optimal mix of theory and experiment, at least as far as the EOS of N*M 
is concerned. 

Our functionals are equally well adapted to the calculation of the EOS of the neutron- 
star crust. For the outer part of the crust, consisting of just electrons and bound nuclei, 
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FIG. 3: (Color online.) Calculated pressures in symmetric nuclear matter (SNM) for functionals 
of this paper (BSkl9 and BSk20 yield almost undistinguishable results from BSk21 on this figure) 
; the box summarizes the range of the measurements reported in Ref. [ijj]. 



we used the HFB-19, HFB-20 and HFB-21 mass models 



161 ]. whenever experimental masses 



are not available. These are based entirely on the BSkl9, BSk20 and BSk21 parameters 



respectively, and were constructed to run from one drip line to the other [3j. Given not 
only the precision fit that these models give to the measured masses but also the constraints 
to NeuM imposed on the underlying forces, no other published mass model will make more 
reliable predictions for the experimentally inaccessible neutron-rich nuclei that dominate the 
outer crust. (For the same reason these mass models are well adapted to the study of the 



r-process of nucleosynthesis 17], ]l8|.) For the inner crust, which consists of neutron-proton 



clusters embedded in a neutron vapor and electron gas, we have likewise calculated the EOS 
with each of the BSkl9, BSk20 and BSk21 functionals, using the ETFSI (extended Thomas- 
Fermi plus Strutinsky integral) method 19j. Our confidence in the use of our functionals in 
this region derives not only from the NeuM constraints to which they have been subjected 
but also from the precision fit to masses, which means that the presence of inhomogeneities 
and of protons is well represented. 

For each of these three Skyrme functionals we thus have at our disposal the zero- 
temperature EOS for all parts of neutron stars, and hence can solve the Tolman- 
Oppenheimer-Volkoff (TOV) equations describing the global structure of spherical non- 
rotating neutron stars 20|, |21|. Solving these equations for different values of the bary- 
onic central density, n cen , gives us the total gravitational mass Ai as a function of the 
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TABLE I: For each functional are shown i) limiting baryonic density for causality, ii) maximum 
neutron-star mass (the quantities in parentheses refer to the corresponding realistic forces) hi) 
corresponding radius and iv) corresponding central baryonic density. (Neutron-star results for 



SLy4 taken from Ref. 



23] 



Force 


ricaus (fm -3 ) 


Mmax/M® 


R (km) 


Ucen (fuT 3 ) 


BSkl9 


1.45 


1.86 (1.84 [7]) 


9.13 


1.45 


BSk20 


0.98 


2.15 (2.20 [9]) 


10.6 


0.98 


BSk21 


0.99 


2.28 (2.3 [8]) 


11.0 


0.98 


SLy4 


1.23 


2.05 


9.99 


1.21 



circumferential radius R. Numerical calculations have been carried out using the LORENE 
library |22j. The maximum possible value of Ai, A4 max , along with the corresponding radius 
R and central density n cen , are given in the last three columns of Table [J for each of our 
three functionals. The second column of this table shows n caus , the maximum value of the 
baryonic density in N*M for which causality holds, i.e., for which the velocity of sound is 
lower than the velocity of light, c (see Section VB of Ref. 3] for details). It will be seen 
that for functionals BSkl9 and BSk20 the central density n cen for maximum mass is at the 
limit imposed by causality; in fact if we had waived this constraint we would have found 
slightly higher (and unphysical) values of A4 max and n cen (the maximum mass would have 
been essentially the same for BSkl9, whereas it would have increased to M. m ax = 2.17 \M Q 
for BSk20). On the other hand, in the case of functional BSk21, which has the stiffest EOS, 
the requirement of causality is satisfied for all possible masses. Table [I] also shows (in paren- 
theses), for each of BSkl9, BSk20 and BSk21, the value of M. max that has already been 
published for the corresponding realistic calculation, i.e., WFF, APR and LS2, respectively. 
The close agreement between corresponding values of Ai max is a measure of how well we 
have reproduced the realistic forces with our Skyrme functionals. 

Comparing the values of M. m ax shown in Tabled with the recently measured value of 1.97 

n 

± 0.04 M G for the mass of pulsar PSR J1614-2230 [24] shows that the EOS of BSkl9 for 
N*M and NeuM is definitely too soft. Even taking account of the fact that t his p ulsar rotates 



with a frequency of 317 Hertz, we have found using the LORENE library 
stationary configurations for rigidly rotating neutron stars, see, e.g. Ref. 



(considering 
25() that M max 
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FIG. 4: (Color online.) Proton fraction Y p in neutron-star matter for functionals of this paper. The 
dashed line indicates the threshold value of Y p for the direct Urea process to be allowed. Dotted 
curves denote supraluminal EOS. 



is raised only by 0.005 Mq, which is insufficient to change this conclusion. On the other 
hand, functionals BSk20 and BSk21, along with SLy4, all provide EOSs that are sufficiently 
stiff to support PSR J1614-2230. 

Referring to Fig. HI we see that BSk21 is the only one of our functionals for which the 
proton fraction can exceed the critical threshold for a direct Urea process to occur, and then 
only for baryon density n > 0.45 fm -3 , i.e., only in neutron stars whose gravitational mass 
M. > 1.59A4 & . As a consequence, our three functionals are compatible with the constraint of 
Ref. 26] that no direct Urea process should occur in neutron stars with masses in the range 
1 - 1.5 Ai & . However, we now see that recently accumulated nuclear-mass data strongly 
favor BSk21. 

Nuclear Masses. Since we fitted our forces to the data of the 2003 AME many new mass 
data have been accumulated. In particular, a new unpublished AME became available in 
2011 23]- This contains 154 new measurements in the interval N and Z > 8, while 9 masses 
that were previously listed as "measured" have now been demoted to an "estimated" status; 
in accordance with our previous practice we do not consider such nuclei. 

The first three lines of Table HT1 show the rms and mean deviations of our HFB-19, HFB-20 
and HFB-21 mass models with respect to all measured masses in the interval N and Z > 8 
in the 2003 AME, while the next three lines show the corresponding results for the subset 
consisting of nuclei with a neutron separation energy of less than 5.0 MeV. A slight favoring 
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of HFB-21 over the other two models is already apparent in this neutron-rich domain, the 
region of interest in neutron-star applications. This trend is strengthened somewhat in the 
next six lines, where we show the corresponding results for the 2011 AME. Accordingly, 
in the next three lines (13 to 15) we look at the deviations with respect to just the 154 
new data points: the discrimination in favor of HFB-21 is now quite marked, although both 
neutron- rich and proton-rich nuclei are included in this set. Of these 154 new nuclei 73 lie 
on the neutron-rich side of the stability line, and we show the corresponding deviations in 
lines 16 to 18 of Table ITTl There is now an even stronger selection of HFB-21. 

Not included in either AME are the 41 unpublished mass measurements of highly neutron- 
rich nuclei by Matos 28). The 2011 AME lists 24 of these nuclei, and although the Matos 
data were not used the agreement is very good. We consider the subset of the remaining 17 
Matos data in the last three lines of Table [Til where we see further confirmation of a strong 
preference for HFB-21. 

Conclusions. Our three nuclear mass models HFB-19, HFB-20 and HFB-21, which are 
based on the functionals BSkl9, BSk20 and BSk21, respectively, give equally good fits to 
the nuclear mass data appearing in the 2003 AME 12]. However, since 2003, and since our 
mass models were fitted, mass measurements have begun to be made sufficiently far from 
the stability line to discriminate between our models. Specifically, model BSk21 is strongly 
favored by the new data. 

To complement this conclusion we have solved the general-relativistic equations of rotat- 
ing and non-rotating stars for our three functionals and find that BSkl9 is too soft at high 
densities to support neutron stars as heavy as the heaviest yet observed, while BSk20 and 
BSk21 will. This corroboration between measurements of nuclear masses and measurements 
of neutron-star masses is gratifying, but it should not be concluded that nuclear-mass mea- 
surements are now tying down the high-density EOS of N*M, since with a more elaborate 
Skyrme parametrization (e.g., by taking several forms of each of the £3,^4 and £5 terms) it 
might be possible to have an equally good fit to all the mass data while making different 
extrapolations to high-density N*M. Nevertheless, at the very least we recognize that in 
BSk21 we have a functional that performs well over a wide range of densities, giving a preci- 
sion fit to nuclear masses, and being compatible with what is now known about neutron-star 
masses. 

We would like to thank M. Matos for sending us his thesis results, and G. Audi, P. Haensel 
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TABLE II: Rms and mean deviations between mass models and different data sets (see text). 



Data set 


Number of nuclei Mass model 


rms (MeV) 


e (MeV) 


2003 AME 


2149 


HFB-19 


0.583 


-0.038 


« 




HFB-20 


0.582 


0.021 






HFB-21 


0.577 


-0.054 


2003 AME n-rich 


185 


HFB-19 


0.803 


0.243 






HFB-20 


0.790 


0.217 






HFB-21 


0.762 


-0.086 


2011 AME 


2294 


HFB-19 


0.595 


-0.010 






HFB-20 


0.593 


0.051 






HFB-21 


0.574 


-0.031 


2011 AME n-rich 


224 


HFB-19 


0.799 


0.261 






HFB-20 


0.788 


0.255 






HFB-21 


0.730 


-0.060 


2011 new 


154 


HFB-19 


0.824 


0.310 






HFB-20 


0.803 


0.375 




)) 


HFB-21 


0.681 


0.185 


2011 new n-rich 


73 


HFB-19 


0.948 


0.391 






HFB-20 


0.893 


0.408 






HFB-21 


0.735 


0.063 


2007 Matos 


17 


HFB-19 


0.912 


0.757 






HFB-20 


0.878 


0.713 




•!1 


HFB-21 


0.536 


0.216 



and D. Lunney for valuable communications. The financial support of the FNRS (Belgium), 
the NSERC (Canada) and CompStar (a Research Networking Programme of the European 
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